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Abstract 

We prove necessary optimality conditions, in the class of continuous functions, for varia- 
tional problems defined with Jumarie's modified Riemann-Liouville derivative. The fractional 
basic problem of the calculus of variations with free boundary conditions is considered, as well 
as problems with isoperimetric and holonomic constraints. 

Keywords: fractional calculus, Jumarie's modified Riemann-Liouville derivative, nat- 
ural boundary conditions, isoperimetric problems, holonomic constraints. 

1 Introduction 

There exists a vast literature on different definitions of fractional derivatives. The most popular 
ones are the Riemann-Liouville and the Caputo derivatives. Each fractional derivative presents 
some advantages and disadvantages (see, e.g., [HSIIMIHO] ) ■ The Riemann-Liouville derivative of a 
constant is not zero while Caputo's derivative of a constant is zero but demands higher conditions 
of regularity for differentiability: to compute the fractional derivative of a function in the Caputo 
sense, we must first calculate its derivative. Caputo derivatives are defined only for differentiable 
functions while functions that have no first order derivative might have fractional derivatives of 
all orders less than one in the Riemann-Liouville sense [55] . 

Recently, Guy Jumarie (see P^H?T][2^26II28 '.'29' ) proposed a simple alternative definition to the 
Riemann-Liouville derivative. His modified Riemann-Liouville derivative has the advantages of 
both the standard Riemann-Liouville and Caputo fractional derivatives: it is defined for arbitrary 
continuous (nondifferentiable) functions and the fractional derivative of a constant is equal to zero. 
Here we show that Jumarie's derivative is more advantageous for a general theory of the calculus 
of variations. 

The fractional calculus of variations is a recent research area much in progress. It is being 
mainly developed for Riemann-Liouville (see, e.g., [JEIEIIIIIIIZIISI] ) and Caputo derivatives (see, 
e.g., [2l[8l[Tni[18l[34l[36|). For more on the calculus of variations, in terms of other fractional 
derivatives, we refer the reader to [1^11151150] and references therein. 

As pointed out in the fractional calculus of variations in Riemann-Liouville sense, as it 
is known, has some problems, and results should be used with care. Indeed, in order for the 
Riemann-Liouville derivatives aD"y{x) and xD'^y{x) to be continuous on a closed interval [a, 6], 
the boundary conditions y[a) — Q and y{b) — Q must be satisfied ^39 . This is very restrictive when 
working with variational problems of minimizing or maximizing functionals subject to arbitrarily 
given boundary conditions, as often done in the calculus of variations (see Proposition 1 and 
Remark 2 of 7 ). With Jumarie's fractional derivative this situation does not occur, and one can 
consider general boundary conditions y{a) — ya and y{b) — j/f,. 
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The paper is organized as follows. In Section [5] we state the assumptions, notations, and the 
results of the literature needed in the sequel. Section|3]reviews Jumarie's fractional Euler-Lagrange 
equations 25 . Our contribution is then given in Section 21 in H4.ll we consider the case when no 
boundary conditions are imposed on the problem, and we prove associated transversality (natural 
boundary) conditions; optimization with constraints (integral or not) are studied in sections i )4.2l 
and §4.31 Finally, in Section [5] we explain the novelties of our results with respect to previous 
results in the literature. 



2 Preliminaries on Jumarie's Riemann— Liouville derivative 

Throughout the text / : [0, 1] — > M is a continuous function and a a real number on the interval 
(0, 1). Jumarie's modified Riemann-Liouville fractional derivative is defined by 

/^"'(.) = ^|:f(.-r"(/W-/(0))ci.. 

If /(O) = 0, then f^"^ is equal to the Riemann-Liouville fractional derivative of / of order a. We 
remark that the fractional derivative of a constant is zero, as desired. Moreover, /(O) = is no 
longer a necessary condition for the fractional derivative of / to be continuous on [0, 1]. 
The {dt)°' integral of / is defined as follows: 

f{t){dtr = a I {x-tr-\f{t)dt. 

Jo 

For a motivation of this definition, we refer to |19) . 

Remark 2.1. This type of fractional derivative and integral has found applications in some phys- 
ical phenomena. The definition of the fractional derivative via difference reads 

f(»){x)=\im^^^P^, 0<a<l, 
■' ^ ' hio h"" 

and obviously this contributes some questions on the sign ofh, as it is emphasized by the fractional 
Rolle's formula f{x + h) = f{x) + h°'f^°'')(x). In a first approach, in a realm of physics, when h 
denotes time, then this feature could picture the irreversibility of time. The fractional derivative 
is quite suitable to describe dynamics evolving in space which exhibit coarse-grained phenomenon. 
When the point in this space is not infinitely thin but rather a thickness, then it would be better 
to replace dx by (dx)"' , < a < 1, where a characterizes the grade of the phenomenon. The 
fractal feature of the space is transported on time, and so both space and time are fractal. Thus, 
the increment of time of the dynamics of the system is not dx but (dx)"' . For more on the subject 
see, e.g., !M[MiM\E3^ ■ 

Our results make use of the formula of integration by parts for the (dx)" integral. This formula 
follows from the fractional Leibniz rule and the fractional Barrow's formula. 

Theorem 2.2 (Fractional Leibniz rule |22j). // / and g are two continuous functions on [0,1] , 
then 

(/(x)5(x))(") = (/(x))(").g(x) + f{x){9{x)Y'^\ (1) 
Kolwankar obtained the same formula ([ij by using an approach on Cantor space |31j . 
Theorem 2.3 (Fractional Barrow's formula 24 ). For a continuous function f, we have 



/(")(t)(dt)(")=a!(/(x)-/(0)). 



where al — T{1 + a). 
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From Theorems 12.21 and 12.31 we deduce the fohowing formula of integration by parts: 



u^°''>{x)v{x){dx)" {u{x)v{x)y°'^dx)°' ^ u{x)v^°'\x) {dx)" 

Jo ^ Jo 

^ al[u{x)v{x)]l - [ u{x)v^°'\x) {dx)°' . 



It has been proved that the fractional Taylor series holds for nondifferentiable frmctions. See, 
for instance, [53]. Another approach is to check that this formula holds for the Mittag-Leffler 
function, and then to consider functions which can be approximated by the former. The first term 
of this series is the Rolle's fractional formula which has been obtained by Kolwankar and Jumarie 
and provides the equality d°'x{t) = a\dx{t). 

It is a simple exercise to verify that the fundamental lemma of the calculus of variations is 
valid for the (dx)" integral (see, e.g., [U for a standard proof): 

Lemma 2.4. Let g he a continuous Junction and assume that 

r-l 

g{x)h{x) (dx)" = 
for every continuous function h satisfying h{0) = h{l) — 0. Then g = 0. 



3 Jumarie's Euler— Lagrange equations 

Consider functionals 



J{y)= / L(x,y{x),y^"\x)) {dx)" 



(2) 

defined on the set of continuous curves y : [0,1] — > M, where L(-,-,-) has continuous partial 
derivatives with respect to the second and third variable. Jumarie has addressed in |25| the basic 
problem of calculus of variations: to minimize (or maximize) ^7, when restricted to the class of 
continuous curves satisfying prescribed boundary conditions ?;(0) = j/o smd j/(l) — yi. Let us 
denote this problem by (P). A necessary condition for problem (P) is given by the next result. 

Theorem 3.1 (The Jumarie fractional Euler-Lagrange equation [25] — scalar case). If y is a 
solution to the basic fractional problem of the calculus of variations (P), then 



dL dL 



dy dx" 



= 



(3) 



is satisfied along y, for all x G [0, 1]. 

Definition 3.2. A curve that satisfies equation (|3]) for all x G [0, 1] is said to be an extremal for 
J. 

Example 3.3. Consider the following problem: 



J{y) 



_r(a + 1) 
subject to the boundary conditions 

2/(0) = 1 and y{l) = 2 
The Euler-Lagrange equation associated to this problem is 



{dxT 



extremize 



r(a + l) 



2x° 



r(a + 1) 



(a) _ 2a;" 



0. 



(4) 



Let y = x"' + 1. Since y*-"^ = T{a + 1), it follows that y is a solution of (|4]). We remark that the 
extremal curve is not differentiable in [0, 1]. 
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Remark 3.4. When a = 1, we obtain the variational functional 

J{y)= f\x{yr-2xyfdx. 
Jo 

It is easy to verify that y = x + I satisfies the (standard) Euler-Lagrange equation. 

We now present the Euler-Lagrange equation for functionals containing several dependent 
variables. 

Theorem 3.5 (The Jumarie fractional Euler-Lagrange equation PS] — vector case). Consider 
a functional J , defined on the set of curves satisfying the boundary conditions y(0) = and 
2/(1) = 2/1, of the form ^ 

J(y)= I L[x,y{x),y^'^Hx)) [dxT , 

where y ~ (j/i, . . . , t/„), y'"-* — {y^\ . . . , yr"^), and yk, k ~ 1, . . . , n, are continuous real valued 
functions defined on [0,1]. Let y be an extremizer of J. Then, 

^ (.,.(.),,(")(.)) --g,^ (.,.(.),,(")(.)) =0, 
k = 1, . . . , n, for all x 6 [0, 1]. 



4 Main results 

We give new necessary optimality conditions for: (i) functionals of form ^ with free boundary 
conditions (Theorem 14. ip : (ii) fractional isoperimetric problems of Jumarie (Theorem 14.31) : and 
(iii) Jumarie- type problems with subsidiary holonomic constraints (Theorems 14.41 and 14. 5p . 



4.1 Natural boundary conditions 

The problem is stated as follows. Given a functional 

r-l 



J{y)= / L(x,y{x),y^^Hx)) [dxT , 



where the Lagrangian L{-, •, •) has continuous partial derivatives with respect to the second and 
third variables, determine continuous curves y : [0, 1] — !■ M such that J has an cxtrcmum at y. 
Note that no boundary conditions are now imposed. 



Theorem 4.1. Let y he an extremizer for J . Then y satisfies the Euler-Lagrange equation 

dL dL ^ 



dy dx" 

for all x e [0,1] and the natural boundary conditions 



dL 



= and ^, , 



0. 



(5) 



(6) 



Proof. Let h be any continuous curve and let j{e) = J{y + eh). It follows that 
= 



dL d°' dL 



dy dx"' dy(°''> 



h{x) (dx)" + a! 



dL 



dyi^) 



h{x) 



x=0 



4 



If we choose curves such that h{Q) — h{l) — 0, we deduce by Lemma 12.41 the Euler-Lagrange 
equation ([5]). Then condition 

dL 



_dy(°' 



■h{x) 



x=l 







x=0 



must be verified. Picking curves such that h{{)) = and h{l) ^ 0, and others such that h{l) = 
and h{Q) ^ 0, we deduce the natural boundary conditions □ 

If one of the endpoints is specified, say y(0) = y^, then the necessary conditions become 

dL (P dL 



and 



dL 







dy 



(a) 



= 0. 



x = l 



Example 4.2. Let J he given by the expression 



The Euler-Lagrange equation associated to this problem is 



and the natural boundary conditions are 



y(-Hx) 



y/l + y(°'\xy 



= and 



= 



y/l + yi°'){xf 



= 0. 



x=0 



Since y^"^ = if y is a constant function, we have that any constant curve is a solution to this 
problem. 



4.2 The isoperimetric problem 

The study of isoperimetric problems is an important area inside the calculus of variations. One 
wants to find the extremizers of a given functional, when restricted to a prescribed integral con- 
straint. Problems of this type have found many applications in differential geometry, discrete and 
convex geometry, probability, Banach space theory, and multiobjective optimization (see [4l[5ll33) 
and references therein). We introduce the isoperimetric fractional problem as follows: to maximize 
or minimize the functional 

J{y)^ f\{x,y{x),y^"Hx)){dxr 
Jo 

when restricted to the conditions 

g{y)^ f' f(x,y(x),y^''Hx)){dxr ^K, K eR, (7) 
Jo 

and 

y(0)-2;oandy(l) = yi. (8) 

Similarly as before, we assume that L(-, •, •) and /(•, •, •) have continuous partial derivatives with 
respect to the second and third variables. 
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Theorem 4.3. Let y be an extremizer o] J restricted to the set of curves that satisfy conditions 
([7]) and ([5]). If y is not an extremal for Q. then there exists a constant A such that the curve y 
satisfies the equation 

dF dF 
dy dx" ~ 

for all X G [0, 1], where F = L — A/. 

Proof. Consider a variation curve of y with two parameters, say y{x) + eihi{x) + £2/12(2;), where 
hi and /12 are two continuous curves satisfying hi{0) = hi{l) — 0, i — 1,2. Consider now two new 
functions defined in an open neighborhood of zero; 

j(ei, £2) = J{y + eihi + £2/12) and ^(£1, £2) = G{y + eihi + £2/12) - K. 
Thus, g(0, 0) = and 



dg 

de2 



(0,0) 



a/ d" df 
dy dx" 



h2{x) {dxT. 



Since y is not an extremal for Q, there must exist some curve ft.2 such that dg/de2{0, 0) ^ 0. By the 
imphcit function theorem we may write g{ei, £2(£i)) — for some function £2(-) defined in an open 
neighborhood of zero. Because (0, 0) is an extremum of j subject to the constraint g{ei, £2) = 0, 
and V(7(0,0) ^ 0, by the Lagrange muhipher rule there exists some real A verifying the equation 



V(j-A5)(0,0)=0. 



In particular. 



= 



dei 



(0,0) 



dL d" dL 



A 



df d" df 



dy rfx" V 9y dx"" i9y(") 

By the arbitrariness oi hi and by Lemma |2.4[ we have 

aL d" dL ^/df_ df 

dy dx"" dy^"'^ \dy dx"" dy'-"^ 

Introducing F = L — Xf we deduce that y satisfies the equation 

dF d°' dF 
dy dx"' dy^"') 

for aU x e [0,1]. 



hi{x) (dx)". 



0. 







□ 



4.3 Holonomic constraints 



In i )4.2l the subsidiary conditions that the functions must satisfy are given by integral functionals. 
We now consider a different type of problem; find functions yi and 2/2 for which the functional 



J{yi,y2)^ J L(x,yi{x),y2{x),y[°'\x),yi°'\x)^ (dx)" 

has an extremum, where the admissible functions satisfy the boundary conditions 

(yi(0), 2/2(0)) = (y°,y2") and (2/1(1), ^2(1)) - {yly^), 
and the subsidiary holonomic condition 

g{x,yi{x),y2(x)) = 0. 



(9) 



(10) 



(11) 
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Theorem 4.4. Given a functional J as in defined on the set of curves that satisfy the 

boundary conditions (1101) and lie on the surface (1111) . let (j/1,2/2) be an extremizer for J . If 
dg/dy2 7^ for all x G [0, 1] , then there exists a continuous function X{x) such that 1/2) satisfy 
the Euler-Lagrange equations 

dF -0 (12) 

" d^ Qyi^) ^ ' 

k — 1,2, for all x G [0, 1], where F — L — Xg. 

Proof. Let y = (yi, 2/2), e be a real, and 2/2(2:)) = y{x) + eh{x) with h{x) = ihi{x), h2{x)) 

a continuous curve such that /i(0) — h{l) — (0, 0). By the impHcit function theorem it is possible 
to solve the equation g(x,yi{x),y2{x)) = for h2, i.e., to write /12 = h2{e,hi). Let j(e) — 
J{yi{x),y2(x)). Then /(O) = and so 



|i..(.) + ^.l")(.) + ^.2(.) + ^4"^(.)j idx) . 

(13) 



On the other hand, we required that {yi{x),y2{x)) satisfy the condition (|lip . This means that 
g{x,yi{x),y2{x)) — 0, and so it follows that 

= ■^gix,yi{x),y2ix)) 



e=0 



-^hi{x) + -;—h2(x) 

dyi dy2 



Thus, we may write /12 in the form 



dg_ 

dyi 

dg_ 

dv2 



h2{x) = -^hi{x). 



Define A as follows: 



dL d° dL 



\{x) = j^^- (14) 

dy2 

Then, we can rewrite equation (jl3p as 

= - -^^4- - X{x)p-\ h,{x) (dxr. 

Jo \dyi dx^dy[^^ ^ ' dyi j ' ^ ' 

By Lemma |2.4[ and since hi is an arbitrary curve, we deduce that 

9L d'' dL Bg 

- X(x)Tr- = 0- (15 

dyi dx^ dy'f^ ^ 'dyi ' 

Let F = L — Xg. Combining equations ((Ti)) and (fT5)) we obtain formula (IT^ . This completes the 
proof of the theorem. □ 

We now state (without proof) our previous result in its general form. 
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Theorem 4.5. Let J be given by 

J{y)= f\{x,y{x),y^'^\x)){dxr, 
Jo 

where y — and y^"^ — . . . , yi"''), such that yk, k = l,...,n, are continuous 

functions defined on the set of curves that satisfy the boundary conditions y(0) — yo and y{\) = yi 
and satisfy the constraint g{x, y) = 0. If y is an extremizer for J , and if dg/dyn ^ for all 
X G [0,1], then there exists a continuous function \(x) such that y satisfy the Euler-Lagrange 
equations 

OF d"' OF _ 
for all x e [0, 1], where F = L — Xg. 



5 Comparison with previous results in the hterature 

The fractional variational calculus, dealing with Jumarie's modified Riemann-Liouville derivative, 
is still at the very beginning. Results available in the literature reduce to those in [3S] (published 
in 2009) and [3ll32l (published in 2010). In {25" the basic problem of the calculus of variations with 
Jumarie's Riemann-Liouville derivative is stated for the first time. The fractional Euler-Lagrange 
equation is obtained under the presence of external forces, and the implication of Jumarie's frac- 
tional calculus in physical situations discussed in detail, in particular Lagrangian mechanics of 
fractional order. A fractional theory of the calculus of variations for multiple integrals, in the 
sense of Jumarie, is proposed in 3 : fractional versions of the theorems of Green and Gauss, 
multi-time fractional Euler-Lagrange equations, and fractional natural boundary conditions are 
proved. As an application, the fractional equation of motion of a vibrating string is investigated [3]. 
In [32] , Euler-Lagrange necessary optimality conditions for fractional problems of the calculus of 
variations which are given by a composition of functionals, within Jumarie's fractional calculus, are 
proved. Optimality conditions for the product and the quotient of Jumarie's fractional variational 
functionals are obtained as particular cases |32) . 

In the present paper we develop further the theory initiated by Jumarie in 2009 [25' , by consid- 
ering other fundamental problems of variational type than those treated before in [3,25,32]. First 
we consider the general form of optimality, when we do not impose constraints on the boundaries, 
i.e., when y{0) and/or y(l) are free (Section l4.ip . If y(0) is not specified, then the transversality 

condition ^ . . (o,y{0),y'-"\0)j = complements Jumarie's Euler-Lagrange equation; if y{l) is 

not specified, then the transversality condition ^ (^l,y{l),y^"\l)^ = holds together with 

Jumarie's Euler-Lagrange equation. We then study fractional variational problems via Jumarie's 
modified Riemann-Liouville derivative under the presence of certain constraints, thus limiting the 
space of continuous functions in which we search the extremizers. Two types of problems were 
considered here for the first time in the literature: fractional isoperimetric problems where the 
constraint also contains Jumarie's modified Riemann-Liouville derivative (Section 14. 2[) : and vari- 
ational problems of Jumarie type subject to holonomic constraints (Section 14. 3|) . The results of 
the paper are valid for extremizers which are not differentiable. 

As future work, we plan to study nondifferentiable fractional variational problems under the 
presence of nonholonomic constraints, which depend on Jumarie's modified Riemann-Liouville 
derivative, and possible applications in physics and economics. 
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